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Abstract: It is very easy to write down models in which long-lived particles decaying 
to standard model states are pair-produced via Higgs decays, resulting in the signature 
of approximately back-to-back pairs of displaced narrow hadronic jets and/or lepton jets 
at the LHC. The LHC collaborations have already searched for such signatures with no 
observed excess. This paper describes a Monte Carlo method to reinterpret the searches. 
The method relies on (ideally multidimensional) efficiency tables, thus we implore collab¬ 
orations to include them in any future work. Exclusion regions in mixing-mass parameter 
space are presented which constrain portal models. 
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1 Introduction 

One of the primary goals of the LHC programme is to study the properties of the recently 
discovered ss 125 GeV state, h, in great detail [1, 2]. So far the state is consistent with a 
standard model (SM) Higgs. Still, plenty of room remains for new physics, particularly in 
its decays; Br{h —)■ unobserved) < 20% is still allowed for an otherwise SM Higgs [3, 4]. It 
is therefore sensible to consider spectacular Higgs decay channels for which searches might 
have high sensitivity. Of interest in this paper is one such channel: pair production of 
long-lived states which decay on detector length scales [5] . 

Long-lived neutral states decaying to SM particles arise naturally in models with ap¬ 
proximate symmetries. There are four mechanisms often encountered in the literature: an 
approximate enhanced Poincare symmetry [6], which if exact would completely decouple 
a hidden sector; an approximate discrete symmetry, which if exact would produce a sta¬ 
ble lightest particle; a low-energy accidental symmetry, with decays only proceeding via 
a heavy off-shell mediator, and; small mass-splittings. The first mechanism is associated 
with the usual portal models. 

Testable models in which a SM-like Higgs decays to pairs of these long-lived states are 
easy to build. The generic signature is an approximately back-to-back pair of displaced 
narrow hadronic jets and/or lepton jets [7] (although this is certainly not the only pos¬ 
sibility), for which null searches have already been performed by the LHC collaborations 
[8-18]. The peculiarity of the signature presents two complementary challenges: how do 
collaborations present their results in the most model-independent way possible? and how 
do phenomenologists reinterpret the results in the context of their own models? This paper 
describes a simple Monte Carlo method which has implications for both. The short mes¬ 
sage is the following: if the relevant efficiency tables are published, then phenomenologists 
need only take Monte Carlo events and fold in these efficiencies to reinterpret searches. No 
displaced decays need be simulated since decay probabilities are easily calculated. 
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The outline of this paper is as follows. In Section 2 we write down some simple 
models of interest with an emphasis on portal models. In Section 3 we describe a method 
for recasting displaced searches by way of example, reinterpreting two ATLAS searches. 
Exclusion regions in mixing-mass parameter space are presented which apply to portal 
models; they are explicitly applied to the Higgs and vector portals in Section 4. Section 5 
concludes. 


2 Models 


In this section we write down some models of interest which should serve to illustrate (a) 
how easy it is to write down a natural model with long-lived states, (b) how these states 
can inherit very different couplings to SM particles, and (c) how an extended model can 
introduce non-trivial final states. Thereupon it should be clear that displaced searches 
would benefit from a model-independent approach. 

1. SM plus real singlet scalar: With an extra real singlet scalar field S, the SM 
potential becomes 

V =+X(p^(j)(pU + C(l)US^+ C(f>US - + >^sS^+ 0^3 + , (2.1) 


where 4> is the Higgs field. The limits a, (3,^ —?■ 0 and C; C 0 a-re technically natural. 

Once ([> develops a non-zero vacuum expectation value (VEV), the cubic term ^(l>^ct)S 
induces mixing between the SM Higgs doublet state and S to produce mass eigenstates 



( cos p — sin p 
sin p cos p 



( 2 . 2 ) 


where p ~ sinp <C 1 is a mixing angle and h is to be identified with the observed SM- 
like Higgs. Then s decays to SM particles through mixing, with a rate oc leading to 
displaced decays when p is very small. 

If Pg < 0, then to first order in a, /3, and C, the mixing p and the effective hss 
coupling K, are 


tanp 




^ _ 

— m^l ’ 


K = (2.3) 


where v k, 246 GeV. 

Terms odd in S can be forbidden by demanding some symmetry, for example Z 2 
{S —)• —S) or classical scale invariance. Then s = S is stable and a dark matter candidate 
[19] unless > 0 and an effective cubic term is generated by S acquiring a non-zero VEV 
(5) = Vs- In this case, at tree-level and to leading order in Q, 


tan p 


(vvs _ ^ iris 


K = (V 


-|- 4m^ 
— 2m^ 


(2.4) 


This is the usual Higgs portal model [20, 21]. 
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The phenomenology of GeV-scale scalars is detailed in Ref. [22]. We note the significant 
QCD uncertainties in the widths for 2mT^ < rus < 4 GeV. The decay width of the Higgs to 
light scalars is 


T{h —)> ss) 


327rm/i 



(2.5) 


Thus for fixed Br{h —?• ss) and rUs, k is fully determined.^ 

2. Two-Higgs-doublet model (2HDM) plus singlet scalar: After symmetry breaking 

in the GP-conserving 2HDM, one ends up with five Higgs fields, of which the neutrals are 
written h,H,A (see Ref. [23] for details). In the decoupling limit, a—)- 7 r /2 — /3, h can be 
associated with the SM-like Higgs. To avoid tree-level flavour-changing neutral currents, a 
symmetry is usually demanded to ensure that each of (tt^, couple to only one of the 

Higgs doublets [24, 25]. The four arrangements, known as Type I, Type H, Lepton-specific 
and Flipped, qualitatively alter the couplings of, in particular, H and A to the SM fermions 
(see Table 2 of Ref. [23] ). 

In the 2HDM with an extra singlet, if the singlet develops a non-zero VEV then its 
real (imaginary) part can mix with h/H (A) to form a mass eigenstate s (a) which could 
be very light. Depending on their admixture of 2HDM Higgs’, the fermion coupling type, 
and the values of tan/3 and a, s and a could couple dominantly to up-types, down-types, 
or leptons, as demonstrated in Figures 7-10 of Ref. [26]. Additionally, if these states are 
mostly singlet then their lifetimes can be long. 

3. Dark photon: A dark photon 7 ^ can arise if an extra U{l)d gauge group is appended 
to the SM [27]. The SM Lagrangian can be extended by 




1 

4 


7 ?' jr'uv 


2 cos 9\y 


TP' RMI' 
^ uv ’ 


( 2 . 6 ) 


where (B''’') is the U{l)d (hypercharge) field strength, 9w is the Weinberg angle, and 
e is the kinetic mixing parameter. If the 7 ^ is massive, and rnz, the couplings of 

7 rf to SM particles are photon-like and oc e? The lifetime of such a 7 ,^ is oc e“^, and if e is 
small enough the 7 ^ can be long-lived. 

A mass term for 7 ^ can be generated through a Higgs mechanism [30-32] . The simplest 
way is to introduce a complex scalar field with a potential 




(2.7) 


such that S acquires a VEV, leaving a massive state s and giving mass « gdVg to the 
dark photon, where gu is the dark gauge coupling.^ In this case, the widths of the SM-like 

^Note that the requirement As < in implies > 2 x 10“® (cf^)^ Br(h^ss) 2 4 There 

is no such implication in the Eq. 2.3 case. 

^In a general model this need not be the case, e.g. it is possible for the extra U(l)d field to mass mix 
with the SM Z boson independently of e [28, 29]. 

®This scenario is natural in the limit e, (/ —>■ 0 and € —>■ 0. Note that an effective ( is generated at the 
level ~ gdffF X (loop factor) by the kinetic mixing term. 
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Higgs to 7 (i 7 rf and ss are oc and it is easy to obtain branchings of 0(10%) when 2ms 
and/or are less than mh [26, 33, 34], If 2m^^ < mg then s decays promptly via 

s —> 7 d 7 d- Otherwise, it turns out, for ms,m^^ ~ GeV it is possible that e and take 
values which result in 0(10%) Higgs branchings and either 7 ^, or s, or both long-lived (see 
Refs. [35, 36] for some discussion). 

Other models: Here we comment on other contexts in which SM-like Higgs decays 
to long-lived states have arisen in the literature. We note that many of these models also 
predict other long-lived state production mechanisms which can be considered indepen¬ 
dently. 

Long-lived right-handed neutrinos can be pair produced and decay to fermion trilin- 
ears [37-40]. In R-parity violating supersymmetry, the would-be neutralino dark matter 
candidate can be long-lived, decaying to fermion bi- or trilinears (see e.g. [41-45]). The 
neutralinos may be pair produced directly or in a cascade if the spectrum below mh is 
sufficiently complex. Models of WIMP baryogenesis also generically predict long-lived par¬ 
ticles which decay via fermion trilinears [46]. In more complex models with a light hidden 
sector, such as hidden valley models, cascades are common [7, 47-49], producing large 
multiplicity final states often associated with missing energy. In a hidden sector with a 
conhning gauge group one expects bound states [50] . If the confinement scale is ~ mh then 
the SM-like Higgs could decay to hidden hadron pairs or hidden glueballs which then decay 
back to SM particles via a heavy mediator (or the Higgs portal) on collider length scales 
[5, 51]. Such phenomenology is typical of recent “neutral naturalness” models [52-54]. If 
the confinement scale is <C mh then one expects hidden jets, which could result in a very 
large multiplicity of displaced vertices along with missing energy [55-57]. 

3 Method and results 

In the following section we describe a Monte Carlo method for recasting displaced searches 
by way of example. The important point to be made is that phenomenologists cannot 
reliably simulate the detector response to displaced decays, and are therefore reliant upon 
efficiency information provided by the collaborations.^ The recast examples will serve to 
highlight which efficiency information is of most interest. 

3.1 Displaced hadronic jets 

The ATLAS Collaboration has presented a search for the displaced hadronic decays of 
pair-produced long-lived neutral particles in 20.3 fb“^ of data collected at ^/s = 8 TeV 
[9, 10]. They considered pair production via the parton process gg ^ ^ ^ where 

is a scalar particle and ttv is a hidden valley pseudoscalar. The decay of tTv is dominated 
by bb for mjr^ > 10 GeV (the cc and rr decays are subdominant, see their Table 1). No 
excess was observed, and limits were placed on the branching fraction of 4> as a function 
of m-iTv lifetime. Presently we describe a method to reproduce the result. 

^This point was also made (and a similar recast method was used) in Ref. [46]. 
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Validation samples of 55 —>• /i —?• ss — {bh){bh) events in ^/s = 8 TeV pp collisions were 
generated using Pythia 8.180 [58, 59] with the default tune. We took mt = 126 GeV and 
mg = 10, 25,40 GeV to match the ATLAS benchmarks.^ 

The cuts used in the ATLAS analysis are listed in the auxiliary Table 6 of Ref. [10]. 
We recreate them as follows. The pair produced long-lived particles are required to have® 

Et{si) > 60 GeV, Et{s 2 ) > 40 GeV, (3.1) 

where Et = E sin0 is a proxy for the measured transverse energy of the resulting jet, and 
the s subscript indicates p^-ordering. This is a fine approximation except for the (non- 
hadronic) s ^ tt decays with ~ 10% branching. We also demand Ai?(si,S 2 ) > 0.4 to 
ensure well-separated jets; this makes very little impact on the benchmarks considered by 
ATLAS, but will matter as mg —?■ None of the relevant remaining cuts, such as on 

isolation and on electromagnetic fraction, nor the trigger efficiency can be replicated since 
no public tool exists to simulate the detector response to displaced decays.^ Therefore 
the remainder of the analysis necessarily involves the folding in of efficiencies provided by 
ATLAS. 

The CalRatio trigger [62, 63] was used to search for vr^ decays at or beyond the edge 
of the electromagnetic calorimeter. This trigger selects narrow jets with Et > 35 GeV, 
logiQ{EH/E em) > 1-2, and a lack of tracks in the inner detector. The trigger efficiency is 
given as a function of radial (longitudinal) decay position for decays in the barrel (endcap) 
region corresponding to the pseudorapidity region |r/| < 1.5 (1.5 < \rj\ < 2.5) in Figure 1 of 
Ref. [9]. Based on these plots we take the trigger-l-reconstruction efficiency of the trigger 
jet to be non-zero and constant only between 2.0 to 3.5 m in the barrel and 4.0 to 5.5 m in 
the endcap, with a respective ratio of 0.20/0.06. The reconstruction efficiency for the non¬ 
trigger jet is not given, but we take it similarly. By construction, the following quantity is 
then proportional to the trigger/reconstruction probability for a given s of lifetime cr: 

f 0.20P(s, cr) if in barrel, 

e(s,CT) = < (3.2) 

[0.06P(s,ct) if in endcap, 

where P is the probability that a state x decays between Lmin and Lmax, 

P(..^) = -exp(-l^)+exp(-l|n). (3.3) 

with 7 and /3 the relativistic parameters for x. The timing of the s decay is required to 
satisfy At < 5 ns with respect to a /3 = 1 particle. This corresponds to requiring an 
absolute decay distance 

Labs<Y^1.5m^L2T- (3.4) 

®The accuracy of the Monte Carlo for an s of mass rris = 10 GeV decaying directly to bb is questionable, 
nevertheless it is possible to force Pythia to do the decay, and it appears that this is what was done in the 
ATLAS analysis. 

^Selections were made nsing the MadAnalysisS vl.l.lObeta SampleAnalyzer framework [60]. 

^Though some attempts have been made [49, 61]. 
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Thus after the aforementioned selection cuts, each event is weighted by a factor 


W{ct) = e(si, cr)e(s 2 , cr), 


(3.5) 


where we take 



for each s, where 9 is the polar angle from the beam line. After this, the remainder 
of the cuts used in the ATLAS analysis should be largely independent of cr. As such. 


penultimately, we rescale the events (with a common number for all cr and nis) to fit 


ATLAS results; the factor turns out to be Ri 6 x 19.0 pb x20.3 /Ngim, where Ngim 
is the number of simulated events. We ignore for simplicity the additional mg-dependent 
< 10% effect related to sub-dominant but non-zero s ^ tt branching. 

After these requirements we find good agreement for the nis = 10, 25 GeV samples, 
but we overpredict for rus = 40 GeV. This is because as rUg approaches Eg (Eq. 3.1), the s 
decay products spread out and the narrow jet trigger efficiency decreases. To properly take 
this into account we require information on how the efficiency depends on the bb opening 
angle, or equivalently, in the limit Eg 3> m^, the boost. In the absence of such information, 
and in an attempt to capture the physics involved, we demand the following (admittedly 
crude) bound on the opening angle of the bh pair from the leading s:® 


AR{b, b) < 1.5. 


(3.7) 


This cut has been tuned so that our results for rUg = 40 GeV best agree with those of 
ATLAS. Note that the spatial separation of the bb pair as seen by the hadronic calorimeter 
is smaller than such a large AR would normally suggest, since the pair appears late. 

In Fig. la we compare the number of events predicted by our analysis to those of 
ATLAS assuming 100% Higgs branching. Despite the apparent crudeness of some of our 
assumptions, we observe good agreement. The 95% CL limit of 20 events (inferred from 
the ATLAS plots) can be used to obtain a limit on the Higgs exotic branching fraction as 
a function of cr. To obtain the exclusion for alternative masses, Pythia signal samples 
g '5 — 7 - L —>• ss — 7 - ( 66 ) ( 66 ) of varying rUg were fed through our analysis. In Fig. lb we present 
our results as limits on Br{h —)• ss) x Br{s —?■ hadronic)'^ as a function of nig and cr. This 
is also a good approximation for the limit on Br{h —)• 'ydld) ^ hadronic)'^. Given 

the good match to ATLAS, we are confident that our results for 10 GeV < nig < 40 GeV 
are reliable, and for nig > 40 GeV are at least sensible. In Sec. 4, Fig. lb is reinterpreted 
to bound mixing-mass parameter space for the Higgs and vector portals. 

We were fortunate in this analysis because the topology of interest was essentially 
already considered by ATLAS for three benchmark values of nig. This allowed us to 
demonstrate the not obvious fact that much of the cr dependence is taken into account 

®Another option is to reject leading s with boost lower than jcut, forbidding particles which have the 
potential to produce opening angles ~ arccos [l — 2/7^4 + — l)/(mg 7 ;f^t)], where rrix is the decay 

product mass. 
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Figure 1: (a) Predicted number of events for the displaced hadronic jet analysis assuming 
100% Higgs branching. Upper: the results from ATLAS (dashed) and our results (solid) 
for rris = 10, 25,40 GeV beginning left-to-right. Lower: the ratio of our results to those 
of ATLAS, (b) Excluded parameter region for s assuming fixed Br{h —)• ss) x Br{s —>■ 
hadronic)'^] the contours mark branchings of 30%, 20%, 10%, 5%, 3%. 


simply by reweighting events with easily calculated decay probabilities (Eq. 3.3). In the 
region where 7 ^ ~ 1 we saw that there was an additional effect that had to be consid¬ 
ered, related to the boost-dependence of the trigger/reconstruction efficiency. This could 
have been anticipated, since any momentum dependence was already integrated out of the 
efficiency plots provided by ATLAS. For this reason our analysis as it stands cannot be 
reliably reapplied to another model since the overall efficiency will scale non-trivially with 
the (correlated) px distributions of the two long-lived particles. However it should serve 
as a conservative estimate for models with more boosted (on average) long-lived pairs. In 
the next analysis the px dependence of the efficiencies is provided and taken into account. 

3.2 Displaced lepton jets 

In Ref. [13], the ATLAS Collaboration presented the search for a SM-like Higgs decaying to 
a long-lived pair of 0(100 MeV) dark photons in 20.3 fb“^ of data collected at -y/s = 8 TeV. 
The benchmark process considered was 

gg-^h^ fdJdi {fd2ld){fd2^d), (3.8) 

where the are hidden fermions and 7 ^ is the long-lived dark photon, inspired by 
Falkowski-Ruderman-Volansky-Zupan (FRVZ) models [48, 49]. Each dark photon de¬ 
cays at or beyond the outer edge of the electromagnetic calorimeter to either pp or eej'n'K, 
resulting in a muon jet (Type 0) or a narrow jet (Type 2) respectively.® No excess over 

®Note that electrons in the hadronic calorimeter will resemble a narrow hadronic jet. 
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the SM expectation was observed and limits were placed on the exotic branching fraction 
of the Higgs as a function of 7 ^ lifetime. These limits are clearly model-dependent, and it 
is not clear how the phenomenologist might translate them. We describe a simple method 
below. 

A validation sample of FRVZ events was generated in Pythia 8.180 [58, 59] by chang¬ 
ing the properties of in-built particles. We took (m/j, ~ (125, 5, 2,0.4) GeV 

and Br(-jd —)• Type 0) = 0.45 to match the ATLAS benchmark. 

The selection criteria (cuts) for the ATLAS analysis are detailed in Table 1 of Ref. [13]. 
We recreate the analysis as follows. The dark photons are required to satisfy 


\v{ 7 d)\ < 2.5, 


lA(/>(7d,,7rf2)l > 1, (3.9) 


as a proxy for the jet acceptance. The remainder of the analysis necessarily involves the 
folding in of efficiencies. After each 7 ^ decays, the final states are of Type 0-0, 0 - 2 , or 2-2 
in the obvious way. After selection, each event is weighted by a factor 


_ 19.2 pb X 20.3 fb ^ 

1T(ct) — Preco\ldii CT)Pj.^co\nd2i ^'T)^trig\1dii'yd2)i 

firm. 


(3.10) 


where Precoild, ct) is the reconstruction probability^*^ for a 7 ^ of lifetime cr, and Strigi'^di, 7 ^ 2 ) 
is the trigger efficiency given that the event is reconstructed. Eq. 3.10 assumes that the re¬ 
construction probability for each of the lepton jets can be considered independently. Both 
Preco and £trig depend on the event Type, and will be described presently. 

The reconstruction efficiencies for a 7 ^ with transverse momentum pT decaying at a 
length L, Sreco = £reco{PT,L), are provided in the ATLAS auxiliary Tables 1-4 [13] for 
Type 0 and Type 2 jets decaying in the barrel and endcap regions, as defined in Table 1. 
We assume £reco = 0 outside of those L,r],pT regions, which appears to be stricter (and 
therefore more conservative) than the barrel/endcap regions used in the full analysis. Since 
the 'jd are very boosted, we do not require a timing veto. The reconstruction probability 
for each jet is then 


Precoi'Jdt ct) — ^ ^ P{pidi CT)£reco{p^ ^L), 

L bins 

where P^'jd, ct) is given by Eq. 3.3 with 

if in barrel, 
if in endcap. 


{Pmin^ Lf, 


= 


' j^min j^raax 
^ xy-h\n ^ xy-h'm 

sin 9 ’ sin 6 


(3.11) 


(3.12) 


Eor any event involving a Type 2 jet, ATLAS used the previously described CalRatio 
trigger. ATLAS provides the CalRatio trigger efficiency Scah defined as the fraction of jets 
passing the offline selection which also pass the trigger, separately as a function of pT and p. 
Type 0-0 events are collected by the 3mu6_MSonly trigger [12, 13], which requires at least 
three standalone (not combined with an inner detector track) muons with px > Q GeV. 

^°This term includes additional rejection criteria such as inner detector isolation. 
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Type 0 

Type 2 

Barrel 

Endcap 

Barrel 

Endcap 

14 < Lxy/cm. < 780 
\g\ < 0.9 

10 <pr/GeV < 100 

50 < L^/cm < 1400 
1.2 < \g\ < 2.5 

10 <PT/GeV < 100 

150 < Lxy/cm. < 410 
\g\ < 1.0 

20 <pr/GeV < 100 

350 < Ljcm < 650 
1.5 < \g\ < 2.4 

20 <pT/GeV < 100 


Table 1: Definitions of barrel and endcap regions for £reco{PTi L) as defined by ATLAS. 


The efficiency of this trigger is dominated by the granularity of the muon spectrometer; to 
reconstruct three muons at least one of the dark photons must produce a pair of muons 
which have a discernible opening angle. ATLAS provides the efficiency £ 2 , defined as the 
fraction of 7 ^ passing the offline selection and also producing two distinguishable muons, 
separately as a function of pT and rj. The efficiency for detecting at least one muon is 
quoted as e>i = 0.8 (0.9) in the barrel (endcap) region. For our purposes we converted 
these efficiencies, making the assumption of independence, into functions of two variables, 
£cai = £cai{PT,i]) £2 = £2{pt,'ii)- For each event, the trigger efficiency given event 
reconstruction is taken as 


^trigildu 7 ^ 2 ) 


/ 

e>i( 7 di)e 2 ( 7 d 2 ) + e>i( 7 d 2 )e 2 ( 7 <ii) - e 2 ( 7 di)e 2 ( 7 d 2 ) if Type 0-0, 

^ ecaz(7dType-o) if Type 0-2, 

^ EcaZ iridx ) + EcaZ {id^ ) “ ^cal (7di Ycal {id^ ) if Type 2-2, 

(3.13) 


where, in an obvious notation, 6 ( 7 ^) = £{p '^This is not quite a model-independent 
trigger efficiency, since £cai and £2 are derived from a lepton-jet gun event sample, in 
which the are generated uniformly in {pT,'n)j f’ut it serves as a good approximation 
for our purposes. After weighting by reconstruction probabilities, we find that this trigger 
efficiency for the FRVZ sample rescales the number of events by an approximately global 
number, ~ 0.5 for =0.1 cm and ~ 0.3 for = 100 cm. 

In Figure 2a we compare the number of events predicted by our analysis to those of 
ATLAS as a function of cr assuming 100% Higgs branching. The obtained ~ 330 events 
at cr = 4.7 cm is an underprediction compared to the full simulation results of 600 ± 40, 
most likely due to the stricter barrel/endcap regions employed. For alternative lifetimes we 
cannot come up with a physical explanation that could account for the shape discrepancy 
between our curve and the reweighted result of ATLAS. 

Now we wish to reinterpret the ATLAS analysis for h —)• jdld decays predicted by the 
vector portal model described in Sec. 2. A signal sample gg ^ h ^ 'Jd'Jd for = 400 MeV 
was generated in Pythia and fed through our analysis. Fig. 2a shows the total number of 
events predicted as a function of cr, as well as broken down by event Type. More events 
are predicted than in the FRVZ model, and they peak at a lower cr, since on average the 7 ^ 
are more boosted. The 95% CL upper limit of ~ 120 (~ 30) on the total (total excluding 
Type 2-2) number of events can be inferred from the ATLAS plots. These numbers can be 
used along with Fig. 2a to limit exotic Higgs branching fractions for niry^ = 400 MeV. 
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(a) 



(b) 


Figure 2: (a) Predicted number of events in the lepton jet analysis assuming 100% Higgs 
branching, = 400 MeV, and Br{'ya —^ Type 0) = 0.45. Upper: the FRVZ model 
results from ATLAS (dashed) and our results (solid) for all event Types (blue/upper) and 
excluding Type 2-2 events (red/lower). Also shown are the number of events expected for 
the h —)• 7 rf 7 rf model in solid black; dotted curves beginning left-to-right show the break 
down in terms of Type 0-0, 0-2, 2-2 events respectively. Lower: the ratio of our FRVZ 
model results to those of ATLAS, (b) Excluded parameter region for 7 ^ assuming fixed 
Br{h — )■ 7 rf 7 rf) X Br{'yci — )• The contours mark branchings of 30%, 20%, 10%, 5%, 

3%. 


For alternative masses, since E ^ m, the properties of the simulated 7 ^ will be ap¬ 
proximately the same but for the boost 7 = E/m^^. If the efficiencies do not change 
significantly surrounding = 400 MeV, which according to ATLAS is at least a good 
assumption for 0.25 < < 1.5 [13], then according to Eq. 3.3 the number of 

events plotted as a function of m^^cT remains invariant. Limits for alternative masses 
can be derived from the = 400 MeV results using this observation. In Fig. 2b 
we present an example exclusion plot derived from Fig. 2a in this way: the limit on 
Br{h — )■ 7 rf 7 rf) X Br{'^d /^/^)^ a function of and cr. This is also a good ap¬ 
proximation for the limit on Br{h —>• ss) x Br{s —)■ /i/i)^. In Sec. 4, Fig. 2b (and related 
limits on branchings to the other event Types) are reinterpreted to bound mixing-mass 
parameter space for the Higgs and vector portals; there the exclusion is extended up to the 
TT threshold m.y^ ^ 3.5 GeV. 


4 Portal limits 

In Figs. 3 and 4 we reinterpret the analyses of Sec. 3 for the Higgs and vector portal models. 
The coloured regions mark the exclusion assuming Br{h —)• ss) = 30%, 20%, 10%, 5%, 3%. 
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Figure 3: Exclusion plot for the real singlet scalar (Higgs) portal (see text). Incremental 
shadings mark areas of non-negligible lifetime. The coloured regions (this analysis) mark 
the exclusions assuming Br{h —> ss) = 30%, 20%, 10%, 5%, 3%. 
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Figure 4: Exclusion plot for the vector portal (see text). Incremental shadings mark 
areas of non-negligible lifetime. The coloured regions (this analysis) and dashed lines 
(CMS Ref. [64]) mark the exclusions assuming Br{h —^ 'Jd'yd) = 30%, 20%, 10%, 5%, 3%. 
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For the Higgs portal Fig. 3, we ask the reader to bear in mind that within the shaded 
region 2771,^ < rUg < 4 GeV the branching fractions and lifetimes are known to be uncertain 
by up to an order of magnitude (see Ref. [22]). We adopt the most recent calculation [65] 
below 1.4 GeV and smoothly interpolate to a perturbative calculation [22] above 2 GeV. 
The limits that can be derived from Fig. 2a (such as in Fig. 2b) are independent of these 
uncertainties, so that there is enough information provided to reinterpret our results for 
alternative branching fractions and lifetimes; in all cases there is excluded parameter space 
with Br{h —>■ ss) < 10%. Limit lines (borrowed from Ref. [22]) are from prompt B decays 
at Belle/BaBar (reds) [66-69], the GHARM beam dump (blue) [70, 71], and LEP (greens) 
[72-74]. Note that prompt searches at LHCb will be more sensistive for 1 < rUs/GeV < 5, 
but non-trivial long-lifetime effects must be carefully taken into account [75]. Very recently, 
LHGb presented a search [76] for long-lived light scalars in B^ —)• K*^s —)> (iF^7r^)(/r+/r“) 
decays; the limit shown in magenta is reinterpreted from their Figure 4 assuming the 
branching expression found in Ref. [77]. Also included are the excluded regions from 
a recent BaBar search for long-lived particles (dark grey) [78], reinterpreted from their 
Figure 3 (1 cm < CTg < 100 cm) assuming Br{B —>■ sXg) ~ 5.6 sin^ P (l ~ [79]. 

For the vector portal Fig. 4, we calculated lifetimes and branchings for < 12 GeV 
using the R{s) = a{e~^e~ —)■ hadrons)/(T(e'''e“ —)• /i''“/r“) values collated by the PDG [80]. 
For > 12 GeV we used those values provided in Ref. [36]. The exclusion lines are 
from a BaBar search (red) [81], supernovae cooling (blue) [82], and various beam dump 
experiments (dark grey) [83-87]. Limits from a recent CMS search [64] for prompt h —)• 
decays are reproduced as dashed lines assuming Br{h —)• jd7d) = 30%, 20%, 10%, 5%, 3%, 
bottom-to-top. 


5 Conclusion 

Searching for the displaced decays of long-lived neutral states is a sensitive way to look for 
exotic Higgs physics. Already, searches at the LHC have probed branching fractions at the 
per cent level. 

In Sec. 2 we demonstrated how simple it is to build natural models with displaced 
Higgs decay phenomenology. Those models taken together serve to emphasise two points: 
(a) that long-lived neutral states can inherit a wide variety of couplings to SM particles, 
and (b) in extended models, many possible production mechanisms and final states exist. 
This motivates collaborations to present their results in as model-independent a way as is 
possible, and in Sec. 3 we sought to explore by example how this might be done. 

Unlike for prompt events, no public tool exists to simulate the detector response to 
displaced decays. Hence phenomenologists wanting to reinterpret displaced searches are 
reliant upon efficiencies provided by the collaborations. In general, the reconstruction 
efficiency for a long-lived particle will depend on its mass m, transverse momentum px, 
lab-frame decay length L, pseudorapidity r/, and its decay mode. Ideally the phenomenol- 
ogist would know the reconstruction efficiency as a function of all five of these parameters 
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together/^ and then the simple Monte Carlo method we have described in Sec. 3, requiring 
no simulation of displaced decays, could be used to determine the reconstruction probabil¬ 
ity for any event. In the interests of simplicity, the dependence on rj is likely to be weak 
enough to be split into barrel and endcap regions, and then the reconstruction efficien¬ 
cies could be provided as three-dimensional {m,pT,L) functions for each final state in the 
barrel/endcap. In certain limits the dependence on one of these variables might even be 
removed. For example, in the limit E ^ m, the efficiency dependence on m for hadronic 
jets is expected to be weak. The trigger efficiency, defined as the probability of triggering 
given reconstruction, could subsequently be taken into account in a similar way. 

Although this efficiency table wishlist was not fully realised for either of the ATLAS 
displaced searches considered in Sec. 3, we were still able to demonstrate the principles of 
a simple Monte Carlo method for reinterpretation, and we used it to constrain parameter 
space of interest for portal models (see also Sec. 4). Our hope is that this paper inspires the 
following take-home message regarding displaced searches: if the relevant multidimensional 
efficiency tables are provided, then phenomenologists will be able to reinterpret searches 
in the context of their own models. 
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